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I. INTRODUCTION 



Brane-world models share a number of characteristic properties which make them a particularly interesting 
starting point for early universe cosmology. One common feature which has stimulated most of the work on brane- 
world cosmology to date is the presence of stress energy localised on the brane leading to inhomogeneous 
additional dimensions. In this paper, we will be concerned with another generic property, namely the possibility 
that branes move as part of the cosmological evolution. While there have been some early proposals for moving 
brane cosmological scenarios pffl the problem has only recently received wider attention |0-n3] . 



Here, we will not be so much concerned with cosmological scenarios but rather with constructing and analysing 
fundamental cosmological solutions with moving branes in a well-defined M-theory context. More specifically, we 
will consider heterotic M-theory [l^-|T^ compactified on a Calabi-Yau three-fold. The resulting five-dimensional 
theory consists of an = 1 bulk supergravity theory on the orbifold jZ^ coupled to two four-dimensional = 1 
theories residing on branes which are stuck to the two orbifold planes M-M. In addition, M-theory five-branes 



can be included in the vacuum of the 1 1-dimensional theory [|l6|,^ 23 1. These five-branes are transverse to the 
orbifold direction and wrap complex cycles in the Calabi-Yau three-fold. Therefore, they appear as additional 
three-branes in the effective D = 5 theory. Unlike the three-branes on the orbifold fix planes, these additional 
branes can move in the orbifold direction and it is this new possibility that we wish to analyse in a cosmological 
context. A five-brane may then collide with one of the boundaries in the course of the cosmological evolution. It 



is then converted into an instanton located on this boundary through a small-instanton transition [^4|-g6| . 

Practically, we will be working in the context of the four-dimensional effective theory, obtained by further 
reduction on the domain- wall vacuum of the D = 5 theory ||T^. The simplest version of this D = A, N = I 
theory contains three chiral superfields, namely the dilaton S, the universal T modulus and the modulus Z which 
specifies the position of the five-brane in the orbifold direction. This effective action constitutes a generalisation 
of the standard low-energy action for Eg, x Eg heterotic string to include the five-bane modulus Z . In this paper, 
we will not explicitly include a non-perturbative superpotential for these moduli but rather focus on the simplest 
possibility of freely rolling fields. 

As we will see, the axions in these superfields can be set to zero consistently. Our problem is then to study 
the cosmological evolution of the dilaton (p measuring the Calabi-Yau volume, the real part of the T modulus /3 
measuring the orbifold size, the five-brane position z and the scale factor a of the universe. It is immediately 
clear from the corresponding effective action, that the fields (p and /3 cannot be set to constants as soon as the 
five-brane moves. Therefore, the evolution of z cannot be studied in isolation and wc must include the standard 
moduli (f) and (3 in our analysis as well. 

Let us briefly summarise our main results. We present a general method of analysing systems with moving 
branes which is based on a cosmological moduli space Lagrangian This allows us to find the most general 
solution for our problem and makes it particularly simple to incorporate a number of generalisations. As usual, 
the solutions we find each exist for the negative as well as the positive time branch. Asymptotically, at either end 
of the negative or positive semi-infinite time ranges, our solutions reduce to the standard rolling radii solutions 
involving a, (3 and 4> only. The five-brane is practically at rest in these limits. However, the early and late 
asymptotic regions correspond to two different rolling radii solutions, in general. The transition between these 
two solutions is generated by the five-brane which only moves during an intermediate period. Interestingly, not all 
rolling-radii solutions can be obtained asymptotically but only those which lead to strong coupling. Hence, while 
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there are standard rolling radii solutions which approach weak coupling at one end of the time range, five-brane 
motion precludes this possibility. Another interesting feature of the solutions is that the five-brane only moves 
a finite (coordinate) distance. For each solution, one can therefore adjust whether or not the five-brane collides 
with one of the boundaries. We present an explicit example, where a negative time-branch solution ends with a 
brane collision accompanied by a small-instanton transition. Our solutions and their generalisations to include the 
effect of a potential provide a generalised framework for negative time branch cosmologies. As such they form the 
common ground between the "traditional" dilaton-driven pre-big-bang scenario psf and the related "ekpyrotic" 
universe scenario [0| where the dilaton is replaced by the five-brane modulus. 



To set the scene, we would first like to present the four-dimensional low-energy effective theory we will be using 
and explain its relationship to the underlying heterotic M-theory. Of particular importance for the interpretation 
of our results is the relation to heterotic M-theory in five dimensions, obtained from the 11-dimensional theory 
by compactification on a Calabi-Yau three-fold. This five-dimensional theory provides an explicit realisation of a 
brane-world theory. 

Our starting point is 11-dimcnsional Horava-Witten theory that is 11-dimcnsional supergravity on the 

orbifold IZ2 x A/10 coupled to two 10-dimensional super- Yang-Mills theories each residing on one of the two 
orbifold fixed planes. Upon compactification of this theory on a Calabi-Yau three- fold to five-dimensions [ p^ , pO[ 
one obtains a five-dimensional = 1 gauged supergravity theory on the orbifold /Z2 x M4. This bulk theory 
is coupled to two D = 4, A'' = 1 theories which originate from the ten-dimensional E% gauge theories and are 
localised on three-branes which coincide with the now four-dimensional orbifold fixed planes. 

An important fact for the present paper is that in compactifying from 11 to five dimensions, additional five- 
branes can be included in the vacuum |^,^ . These five-branes are transverse to the orbifold and wrap (holomor- 
phic) two-cycles within the Calabi-Yau three-fold. Hence, they also appear as three-branes in the five-dimensional 
effective theory. Unlike the "boundary" three-branes which are stuck to the orbifold fix points, however, these 
three-branes are free to move in the orbifold direction. Altogether, the five-dimensional theory is then still a 
gauged A^ = 1 supergravity theory on the orbifold /Z2 x M4. It now couples to the two A^ = 1 theories on the 
boundary three-branes as well as to a series of A^ = 1 theories residing on the additional three-branes. 

The vacuum of this five-dimensional theory is a non-trivial BPS domain wall solution with an associated effective 
four-dimensional theory which describes fluctuations around this vacuum state [|9|. It is this four-dimensional 
effective theory, in its most elementary form, which we will use as a starting point for the present paper. First 
of all, for simplicity, we focus on the case of one additional five- (three)-brane in the bulk. The generalisation 
of our results to a more complicated configuration is straightforward and will be briefly discussed in the end. 
For our purpose, we can truncate the four-dimensional theory to three chiral superfields, namely the dilaton 
S, the universal T modulus and the position modulus Z which specifies the position of the additional brane in 
the orbifold direction. The Kahlcr potential which governs the dynamics of these fields has been computed in 
Ref. |2^,^ and is given by 



where 95 is a constant. Let us briefly explain the meaning of (the bosonic components of) these superfields in 
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II. FOUR-DIMENSIONAL EFFECTIVE ACTION 




(1) 



terms of the fields in the underlying five-dimensional theory. To do this we introduce the six real fields {4>,a), 
{P,X) and (zX) through 

S = e'l' + q^z^e^ - 2i{a - q^z'^x) 

T = e^ + 2ix (2) 

The field (j) originates from the five-dimensional dilaton which is part of the universal D = 5 hypermultiplet. It 
measures the size of internal Calabi-Yau three-fold (averaged over the orbifold) such that the Calabi-Yau volume 
is given by we*^ with some fixed reference volume v. The corresponding dilatonic axion a is the other zero mode 
surviving from the universal hypermultiplet. The field (3 is obtained as the zero mode of the (55) component in 
the D = 5 metric and as such measures the size of the orbifold. More precisely, this size is given by wpe^ with 
a fixed reference size np. The associated axion x is related to the vector field in the five-dimensional gravity 
multiplet. Finally, the field z originates from the world-volume of the additional three-brane and specifies the 
position of this brane in the orbifold direction. We have normalised this field to the interval 

ze[0,l] (3) 

where the two endpoints z — and z — 1 correspond to the locations of the boundary branes. The associated 
axion also originates from the additional three-brane and it is related to the self-dual two-form of the underlying 
five-brane. The constant gs in the above Kahler potential can be written as 

95 = 7reo/35 (4) 

where (3^ is the five-brane charge and eg is defined by 

/ K \2/3 2ttp 

^0 -[j^) ^ • (5) 

Here, k is the 11-dimensional Newton constant while v and np have been introduced above. We note that the five- 
brane charge and, hence, the constant q^ needs to be positive if the compactification is to preserve four-dimensional 
supersymmetry, as we have assumed. 

The above Kahler potential is a lowest order expression in the strong-coupling expansion parameter e defined 

by 

e = eoe"-^ . (6) 

This quantity measures the size of the loop corrections to the four-dimensional effective action, or, from a five- 
dimensional viewpoint, the warping of the domain wall vacuum solution. While the S, T part of the Kahler 
potential does not receive any corrections at first order pof in e, such corrections are expected for the Z-dependent 
part. To date, these corrections have not been computed and we will work with the lowest order expression which 
is valid as long as e is sufficiently small. 

When the five-brane collides with one of the boundaries, that is, if 2: — > or z ^ 1, one expects the appearance 
of new massless states which originate from membranes stretched between the five-brane and its Z2 mirror and 
the five-brane and the boundary. The theory then undergoes a small-instanton transition p^-[2^ , where the bulk 
five-brane is converted into an instanton (or a gauge five-brane) on the boundary. The particle content and other 
properties of the = 1 theory on the relevant boundary can change dramatically under this transition. At 
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present, the dynamics are not well enough understood to be able to follow the cosmological evolution through 
this transition without additional ad-hoc assumptions. Here, we will not attempt to improve on this but we will 
merely analyse the conditions under which a brane-coUision occurs. 

Perturbatively, the superpotential for the moduli fields S, T and Z vanishes but non-trivial contributions are, of 
course, expected at the non-perturbative level pl|-|3^ . While there is a long history of analysing the phenomenol- 
ogy of non-perturbative superpotentials for S and T, particularly those induced by gaugino condensation, only 
recently first efforts have been made to include the modulus Z [|ri|,Q^. While this is an interesting direction 
for future research, here we focus on the simplest case of a vanishing superpotential. In our cosmological context, 
this amounts to finding generalisations of rolling-radii solutions to include the position modulus Z. In general, 
rolling-radii solutions constitute a basic and important class of string cosmological solutions and appear to be the 
natural starting point for a study of moving-brane cosmologies. 

The component action of our model can be straightforwardly computed from the Kahler potential (|l|) together 
with the reparametrisation (H). To fix our conventions, we provide the relevant terms in the supergravity action 
which are given by 

1 

Here K^j = d^K/d^^d^^ is the Kahler metric and = (5, T, Z). In terms of higher-dimensional quantities, 
the four-dimensional Newton constant up can be expressed as 

«p = ■ (8) 

It can be seen that all terms in the resulting component action are at least bilinear in the axion fields cr, x C- 
Therefore, these fields can be consistently set to zero and we can work with an action that is further truncated 
to the fields (f), (3 and z. The component action then reads 

1 



(7) 
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ii? + \d^<t>d^^(t> + \d^.Pd^^P + '^ef'-'^d^zd^'z 



(9) 



It is for this simple action that we wish to study cosmological solutions in the following sections. Due to the 
non-trivial kinetic term for z, solutions with exactly constant (/) or /? do not exist as soon as the five-brane 
moves. Therefore, the evolution of all three fields is linked and (except for setting z = const) cannot be truncated 
consistently any further. Recently, this has also been pointed out in Ref. We remark that this statement 

remains true in the presence of a non-perturbative superpotential. 



III. COSMOLOGICAL SOLUTIONS WITH A MOVING FIVE-BRANE 

Let us now look for cosmological solutions to the field equations derived from the action (^) . For simplicity, we 
assume the three-dimensional spatial space to be flat. Our Ansatz then reads 





(10) 


= 


(11) 


a — q;(t) 


(12) 


/3 = (3{t) 


(13) 


z = z{t) 


(14) 
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Later on, we will indicate how to include three-dimensional spatial curvature. The equation of motion for z can 
be immediately integrated once, leading to 



Z = ^,5-3"+'^-^+-^ , (15) 

where the dot denoted the derivative with respect to r. This result may be substituted back into the rest of 
the system to obtain a closed set of equations for a, /3 and 0. These equations can then be integrated in a 
manner closely analogous to the formalism for anti-symmetric tensor fields in cosmology developed in Ref. pTf . 
We will subsequently use this formalism which makes the general structure of the equations and their solutions 
particularly transparent. Also, several generalisations to be discussed later are most easily incorporated in this 
language. 



A. Equations of motion and their solutions 



It turns out that the remaining equations of motion for our system can written in the elegant form 

dr oa 
-Ea'^^Ga' + E-^U = . 



(16) 
(17) 



Here, the scale factors have been arranged in a vector a = (a,/?, 0) which spans the moduli space of our 
cosmological model. The potential U on this moduli space is given by 

(18) 



U 
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where the vector q is given by q = (0,-1,1). Further, the "Einbein" E is defined hy E = g-'^+'i " with the 
dimension vector d = (3, 0, 0) and we have introduced a metric G = diag (—3, |, |:) . Clearly, these field equations 
can be obtained from the Lagrangian 



£ = -Ea''^Ga' - E-^U 
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(19) 



which encodes the dynamics on our moduli space in a compact form. The presence of the potential C/ is a direct 
consequence of the five-brane motion and its five-brane origin is encoded in the characteristic vector q. In more 
complicated situations, for example when spatial curvature is included, the potential U consists of more than one 
exponential. Even then, the Lagrangian ( p^ can sometimes be solved by employing Toda theory. We will return 
to such generalisations later on. 



However, for now, we are dealing with the simple potential (|18|) containing one exponential due to the moving 
five-brane only. This system can always be integrated and the general solution in the comoving gauge, v = 
written in vector notation, is given by 



a = p., In 



T 



(p/ 



I In 



t - tc 
T 



-<5 



Q!0 



(20) 



d 1 



T 



(21) 



where t is the proper time. Let us discuss the various integration constants in this solution. The time-scales 
and T are clearly arbitrary as are the constants d and zq which parameterise the motion of the five-brane. In 
contrast, the initial and final "expansion powers" Pi and p/ each satisfy the two constraints 

p„Gp„ = , p„ • d = 1 , (22) 

where n = i, f, and, moreover, they are mapped into one another by 

2G-iq 



Pf = 



(23) 



< q,q > 

Here, the scalar product < •, • > is defined hy < a,b >= a^G^^h. The power 6 is defined by 

5 = -q • p,; . (24) 

We note that the map ( p3| ) and the conditions (|2|) are consistent with one another in the sense that p/ as given 
by Eq. (|2^) will always satisfy the constraints (|2^) as long as the vector p^ does. As a consequence, the six 
expansion powers contained in p^ and pj depend on one independent degree of freedom. Finally, the additive 
constants ckq are subject to one further condition 

'(Pq5<q,(l>\ 

(25) 



q • ckq = In 

Altogether, we have seven independent integration constants as expected on general grounds. 



t e 



(26) 



In order for the logarithms in Eq. (|2C| ) to be well defined the range of t should be restricted to 

[— oo,io] , (—) branch 
[to,+oo] , (+) branch 

As usual, we, therefore, find a positive and a negative time branch for each solution. One can easily verily that 
the former always starts out in a past curvature singularity while the latter ends in a future curvature singularity. 



(27) 



A useful observation is that each of the above solutions is invariant under the mapping 

Pi ^Pf , Pf ^Pi 1 -d , zo ^ zq + d . 

At the same time, this map implies a sign change S — > —S as can be seen from Eq. ( p3| ) and (|2^). In order to 
avoid having two different sets of integration constants associated to each solution we can, therefore, fix the sign 
of the parameter S. It proves convenient to adopt the convention 

S > (-) branch 
S <0 (+) branch 



(28) 



Finally, for the sake of concreteness, let us present the explicit form of our solutions by inserting the particular 
values for the various general quantities which we have defined. From Eq. ( |20| ) the scale factors can be written 
in component form as 



a = - In 
3 



t-to 



(3 = pp^i In 



In 



T 

t~to 



T 

t-to 



T 



ao 

+ {Pf3J - Pp,i) In 
6 



t-t. 



T 



t-to 



T 
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1 +/3o 



+ 1 + -^0 



(29) 
(30) 

(31) 



while the expression (|2l|) for the five-brane modulus remains unchanged. We see that both expansion powers for 
the scale factor a are given by 1/3, a fact which is expected in the Einstein frame and can easily be deduced from 
the second equation (|2|). The initial and final expansion powers for /3 and (p are less trivial and are subject to 
the first constraint ( p^ ) which explicitly reads 

^Pl,n+Pl,,i = ^ (32) 
for n = i, /. From Eq. (|2^ ) it follows that they are mapped into one another by 

We note that this map is its own inverse, that is P" = 1, which is a simple consequence of time reversal symmetry. 
The power S is explicitly given by 

S = Pfi.i - P<p,i ■ (34) 
Finally, the condition (|2^) takes the explicit form 

cbo-f3o = \n(^) . (35) 



B. Properties and asymptotic behaviour of the solutions 

Let us now discuss the physical properties of our solutions in some detail. For both branches, the system 
behaves in a particularly simple manner at either end of the allowed t ranges. In these limits it is easy to see 
from our solution as presented in equation (|o|) that the 'Hubble parameters' can be written as, 

H = a ~ (36) 

with constant expansion coefficients p satisfying, 

pGp = , d • p = 1 . (37) 

In view of our convention (^8|), the expansion power p in these relations is given by p = p^ at early times (that is, 
at i — )■ — oo in the (— ) branch and t ~* to in the (+) branch) and by p = p/ at late time (that is, at t —^ to in the 
(— ) branch and t ^ cx) in the (+) branch). The above relations characterise the standard rolling radii solutions 
for a, (3 and cj) which have been used as a basis for pre-big-bang cosmology |3^]. They can be obtained as exact 
solutions of the action (^) when the five-brane is taken to be static. Their appearance as asymptotic limits of our 
full solutions can be understood from the five-brane evolution (pl|). This equation shows that, asymptotically, 
the five-brane does not move and effectively drops out of the dynamics thereby leaving us with freely rolling radii 
in these limits. This demonstrates that our solutions show inflationary behaviour at least in some region of the 
negative time branch. 

Thus, we arrive at the following interpretation for our solutions. At early time, the system starts in the rolling 
radii solution characterised by the initial expansion powers Pi while the five-brane is practically at rest. When 
time approaches |t — to| ^ \T\ the five-brane starts to move significantly which leads to an intermediate period 
with a more complicated evolution of the system. Then, after a finite comoving time, in the late asymptotic 
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region, the five-brane comes to a rest and the scale factors evolve according to another rolling radii solution with 
final expansion powers p / . Hence the five-brane generates a transition from one rolling radii solution into another 
one, a process which is formally described by the map (|3^). To illustrate the situation, in Fig. |l|, we have plotted 
the allowed expansion powers for rolling radii solutions in the pp-p^ plane. The arrows indicate the map between 
early and late expansion powers generated by the five-brane. 
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-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

FIG. 1. It is shown how the ellipse of possible expansion powers {pi3,p^) is mapped into itself under the influence of the 
moving five-brane. Points related by the map may be joined with a straight line of slope —3 as indicated by the arrows. 
There are two fixed points under this map on the ellipse given by S — pp ~ pa =0. These points correspond to solutions 
which are related to separating solutions in five dimensions. In these special cases the five-brane modulus z is constant. 



An obvious question to ask is whether all of the rolling radii solutions represented by the ellipse in Fig. |^ are 
available as asymptotic limits at, say, early time. Perhaps surprisingly, the answer is no. As can be seen from 
Eq. (|8|), the expansion powers at early time are subject to the constraint 

^ I > (— ) branch 

Therefore only the solid (dashed) half of the ellipse is available as the early time limit in the negative (positive) 
time branch. This can also be understood from the consistency requirement that the kinetic energy of /3 and (f) 
be much larger than the kinetic energy of the five-brane modulus asymptotically if the system is to approach a 
rolling radii solution. In fact, it can be easily seen that the ratio of these kinetic energies is proportional to the 
strong coupling expansion parameter e defined in Eq. (|^). For our set of solutions, the asymptotic time evolution 
of this parameter is specified by 

e^|l*-^"C (39) 
[ \i~tQf> late 



where we have defined 



Si = S=-q-pi, Sf ^ -di ^ -q- pf . (40) 

We see that, the sign of the e expansion power is inverted as we go from the early to the late asymptotic 
region. If we pick a solution with the "wrong" sign, for example a Si < solution at early times in the minus 
branch, the share of the five-brane kinetic energy increases as t ~ to ^ —oo which is inconsistent with the rolling 
radii nature of the solutions in this limit. Hence such solutions do not exist and we are restricted to the cases 
which respect (^). As a consequence, e always grows asymptotically. This implies growing loop corrections 
asymptotically or, equivalently, an increasing warping in the orbifold direction. Hence, while there are perfectly 
viable rolling radii solutions which become weakly coupled in at least one of the asymptotic regions, the presence 
of a moving five-brane always leads to strong coupling asymptotically. A similar phenomenon has been observed 
when axion dynamics is added to the simple rolling radii set-up | |3^ , |37[ | . 

We note that there two exceptional solutions, indicated by the solid dots in Fig. ^ which are the fix points 
under the map (^). In other words, for those solutions the early and late expansion powers are identical. It is 
clear from Eq. ( p3| ) that those solutions are characterised hy S — or, equivalently, hy — P4, = ±-^. Hence, 
for those solutions the orbifold and the Calabi-Yau space always evolve at the same rate and the size of the 
loop corrections is time- independent. However, it follows from Eq. (^l|) that also the five-brane is static in this 
case, so that these are really ordinary rolling radii solutions. They correspond to the two separating cosmological 
solutions of five-dimensional heterotic M-theory found in Ref. In conclusion, despite these exceptional cases, 
it remains true that the solutions are always driven to strong coupling asymptotically as soon as the five-brane 
moves. 

Our simple effective action (|^) remains valid only as long as loop and higher derivative corrections are sufficiently 
small and this limits the validity of our solutions in the usual way. Another, perhaps more interesting event leading 
to a break-down of the effective action is the collision of the five-brane with one of the boundaries, that is z — > 0, 1. 
In this context, an important observation from Eq. ( pT| ) is, that the five-brane moves a finite coordinate distance 
given by jdj, only. Consequently, whether or not a brane collision takes place, depends on the choices for the 
arbitrary integration constants d and zq. This may be surprising, perhaps, as one could have expected that a 
collision is unavoidable in a cosmological setting and in the absence of any stabilising potentials. Of course, this is 
related to the system always becoming strong-coupling asymptotically which implies that the non-trivial kinetic 
term for five-brane modulus leads to damping terms in the z equation of motion. We also remark that, for each 
solution, the five-brane can move in either direction, as specified by the sign of d. It can, therefore, collide with 
either one of the boundaries. This is, of course, a direct consequence of our simple effective action (^) which is 
invariant under z —z and does not single out a particular direction. 



IV. AN EXPLICIT EXAMPLE 

We would now like to illustrate our general results by an explicit example. For concreteness, we focus on 
the negative-time branch and consider the solutions with an approximately static orbifold at early time. In the 
context of our formalism, these solutions can be singled out by setting the initial expansion power of the orbifold 



to zero, that is pi^p = 0. From Eq. (32) this implies two possibilities for the initial expansion power of the 



dilaton, namely ^ — i"^- Only one of these two expansion powers, pi^^ — leads to a positive value for 



S = Pi, 13 — Pi,4> = and can, therefore, be realised as the early time limit of one of our solutions in the negative 
time branch. Using the map (p3|), the complete set of initial and final expansion powers for our example is then 
given by 

( \ \ 



I 



P» = 



P/ = 



\ 

3 

1 

V3 

V 7f ) 



(41) 



Recall here that the first entry in those vectors specifies the expansion power of the scale factor a (which always 
equals ^) while the second and third entry refer to /3 and 0, respectively. 




FIG. 2. Time-behaviour of P (upper curve) and (j) (lower curve) for the example specified in the text. 



At early times, \t — to\ ^ |T|, the evolution is basically of power-law type with powers p^. As discussed, 
this happens because at early time the five-brane is effectively frozen at z ~ d + 2o and does not contribute a 
substantial amount of kinetic energy. This changes dramatically once we approach the time |i — to| ^ jT'l- In a 
transition period around this time, the brane moves from its original position by a total distance d and ends up 
at z ~ zq. At the same time, this changes the behaviour of the moduli /? and cf) until, at late time \t\ <^ |T|, they 
correspond to another rolling radii solution with powers controlled by p^. Concretely, the orbifold size described 
by f3 turns from being approximately constant at early time to expanding at late time, while the Calabi-Yau size 
controlled by (j) undergoes a transition from expansion to contraction. This behaviour has been illustrated in 
Fig. I 

Recall from Eq. (|3^) that the expansion powers of the strong coupling expansion parameter e in the initial and 
final asymptotic regions are given by (5„ = pp^n — P4>.n, where n — i,f. We have also seen in general that the 
transition induced by the brane is simply inverts the sign of this expansion power, so that df = —Si. For our 
example, we find 

S^ = 2/V3^ Sf = -2/V3 (42) 
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so that e decreases initially and increases after the transition. Consequently, the solution runs into strong coupling 
in both asymptotic regions t — to — oo and t — to which illustrates our general result. 

Of course, the range of validity of our solution is constrained by loop corrections, as well as higher curvature 
corrections, in a way that depends on the various integration constants involved. If the constants d and zq, 
characterising the five-brane motion, are chosen such that the five-brane does not collide with the boundary there 
is no further constraint. As discussed, this can always be arranged. However, we may also choose d sufficiently 
large so that the five-brane does collides with one of the boundaries. Then the system undergoes a small-instanton 
transition and our classical solution is invalidated. This typically happens at the time |t — ^ol — l^^l when the 
five-brane moves significantly. Depending on integration constants this may happen well before loop and higher- 
derivative corrections become relevant. In Fig. ^ we have shown a particular case for our example which leads to 
brane collision. The five-brane is initially located at d -f zq ~ 0.9 and moves a total distance oi d = 1.5 colliding 
with the boundary at z = at the time \t — tol/l?"! — 1- 




FIG. 3. Time-behaviour of the five-brane position modulus z for the example specified in the text. The boundaries are 
located at z = 0,1 and the five-brane colUdes wtth the z — boundary at \t/T\ ~ 1. 

This represents an explicit example of a negative-time branch solution which ends in a small-instanton brane- 
collision. 

V. GENERALISATIONS AND OUTLOOK 

In this paper, we have presented and analysed the simplest class of cosmological solutions with moving five- 
branes in the context of four-dimensional effective heterotic M-theory. Clearly, there are a number of interesting 
generalisations and refinements some of which we would briefly like to discuss below. 

• Including spatial curvature : Consider the case where, in addition to the elements included in the above 
discussion, we wish to include three-dimensional spatial curvature. One of the benefits of our general 
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formalism is that this can be easily incorporated. One can simply continue to work with a moduli space 
Lagrangian of the form ([l9|). The only change is that another exponential term proportional to exp(qc • a.) 
which accounts for the spatial curvature has to be added to the potential U. The characteristic vector 
Qc of this curvature term is different from the corresponding five-brane vector and is explicitly given by 
Qc = (4, 0, 0). Hence, the five-brane and curvature characteristic vectors are orthogonal with respect to our 
scalar product, that is, < q, qc >= 0. Therefore, the system can be identified as an SU{2)'^ Toda model 
which can be integrated. In a different, but formally similar context, this has been explicitly carried out in 
Ref. 

• Add a perfect fluid : One may want to consider the cosmology of our simple model in the additional presence 
of a perfect fluid. As before, this can be done by simply adding an exponential term to the potential U. 
The characteristic vector for a perfect fluid with equation of state p = wp, where w is a constant, is given 
by qfiuid — (3(1 — uj), 0, 0) which is again orthogonal to the five-brane vector. Hence this situation is again 
described by an SU{2)'^ Toda model. 

• More than one five-brane : We may also easily generalise our results to include more than one five brane. 
This simply involves including more potential terms with q = (0,-1,1). This, of course, simply amounts 
to replacing the expression for the normalisation of the potential U in Eq. (|lj) by a generalisation 
which involves summing over all five-branes. Most of the qualitative features of our solutions will then 
be unchanged. In particular, all five-branes 'move at the same time' collectively generating the transition 
between the two asymptotic regions. However, the initial position and the total distance of the motion can 
be different for the various five-branes. For example, some of the five-branes may collide with a boundary 
while the others do not. Of course the five-branes may also collide with one another. 

• Include axions : Another obvious generalisation to attempt would be to try and include non-trivial axion 
dynamics in our solutions. In more simple string cosmology settings without five-branes this has been 
possible due to an SL{2) symmetry of the string effective action [ ^6|j3^ . This symmetry allowed solutions 
with non-trivial axion dynamics to be generated from systems where the axion fields had been set to zero. 
However, the presence of the five-brane modulus Z complicates the situation significantly. The only parts 
of SL{2) which are obviously respected by the Kahler potential (|l]) are the axionic shift symmetries which 
do not lead new solutions. It is not clear, at present, whether the Kahler potential (|l]) has any further 
symmetries which could be used to generate new cosmological solutions from the ones obtained in this 
paper. 

• Small instanton transition : Finally, an understanding of the dynamics of the small instanton transition 
whereby a five brane is emitted from or absorbed onto one of the fix point branes would be extremely 
useful. It would allow us to describe the system as it evolves through the transition which is vital in order 
to construct "complete" cosmological models. 

• Fluctuations : It is known |^,^ that there exist special rolling radii solutions which, in the context of 
pre-big-bang inflation, lead to perturbations consistent with a Harrison-Zeldovich spectrum. It is, therefore, 
expected that such cases with an acceptable fluctuation spectrum also exist in our generalised framework 
with moving five-branes. We will be analysing this in more detail. 

Another possible direction for future work would be a search for five dimensional moving brane solutions of 
various kinds, including those that correspond to the four dimensional ones found here. It is, of course, straight- 
forward to "oxidise" all our four-dimensional solutions to approximate solutions of five-dimensional heterotic 
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M-theory following the procedure described in Ref. However, finding the exact solutions in five dimensions 
in the presence of additional five-branes is a much more difficult task. We have seen that the five-brane is static 
for the two exceptional solutions with 5 — 0. In this sense, they are trivial extensions oi D — A rolling radii 
solutions. Finding their D = b counterpart, however, is less trivial, since the presence of the additional five-brane 
changes the warping of the D = 5 domain wall, an effect which is not directly visible in Z? = 4. However, one can 
verify that the two separating D = 5 cosmological solutions found in Ref. j^] can indeed be generalised to include 
five-branes. Hence, the D = 5 counterparts of our two exceptional solutions can be obtained as exact solutions. 
More generally, we have shown that there exist only three D = 5 solutions, based on the bulk BPS domain wall 
solution of Ref. |19|, which include a single five-brane. These are the BPS domain wall solution with a single 
five-brane itself and the two aforementioned cosmological solutions. Hence, lifting the other solutions we have 
presented to exact D ^ 5 solutions is not straightforward and requires the study of more complicated metrics in 
the bulk such as those foimd in Ref's [pq,p9l. 
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